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Abstract 

Many physical systems like supersymmetric Yang-Mills theories are for- 
mulated as quantum matrix models. We discuss how to apply the Bethe 
ansatz to exactly solve some supersymmetric quantum matrix models in the 
large- limit. Toy models are constructed out of the one-dimensional Hub- 
bard and t-J models as illustrations. 
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More and more physical systems are formulated as quantum matrix mod- 
els; notable examples are quantum cliromodynamics [0, ||, 0, ^, supersym- 
metric Yang-Mills theory D-branes 0, 0, M-theory |Ty, |1T[ and string 

theory |jl2|, |T^. The dimensions of these matrices can be interpreted as the 
numbers of colors or the numbers of D-branes. Taking the limit as the di- 
mensions going to infinity is also of interest; for instance, it is conjectured 
that the large- limit captures the essence of quantum chromodynamics like 
quark confinement and that M-theory in the infinite momentum frame 
is exactly described by the large- limit of DO-brane quantum mechanics 
p!0| . Certainly if ways are developed to solve quantum matrix models in 
the large- iV limit, we will understand better the physics of these interesting 
systems. 

In this Letter, we will concentrate ourselves on matrix models whose 
Hamiltonians preserve the number of partons. These models can be used to 
describe asymptotic spectra of supersymmetric Yang-Mills theories 0. It is 
found that the ubiquitous Bethe ansatz 0] is helpful in solving some of these 
models, and the Yang-Baxter equation can help us determine if the Bethe 
ansatz is applicable. In previous papers |]15|, [1^ , we provided some examples 
of bosonic quantum matrix models which can be solved by the Bethe ansatz. 
Here, we will give first instances of quantum matrix models which contain both 
fermions and bosons and which are still exactly integrable. Previously, work 
was done to solve some integrals over finite chains of matrices with both 
commuting and anti-commuting entries which describe classical integrable 
systems with supersymmetry (see Ref. |ll7| and the citations therein). From 
this perspective, the implication of our discussion here is that certain path 
integrals over matrices corresponding to quantum integrable systems with 
supersymmetry can be exactly evaluated, too. 

Consider a physical system of bosons with degrees of freedom (other 
than colors or labels for D-branes) and fermions with Kp degrees of freedom 
(again other than colors). (In the context of D-branes, A^ is the number of 
transverse dimensions of the D-branes; in the context of gauge theory, A^ 
gives the number of different values of the longitudinal momentum a gluon 
can take.) Let the vacuum be trivial. Let a^^{k) be an annihilation operator 
for a boson if 1 < /c < A^, or a fermion if A^ + 1 < < A^ + Kp- Here 
Hi and are row and column indices running from 1 to (physically they 
can be color indices or labels for D-branes). The corresponding creation 
operator is aj^^^(A;). They satisfy typical (anti)-commutation relations, the 
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most important ones being 



for 1 < fci, A;2 < A^; and 

for Aij + 1 < /ci, A;2 < As + Ap- 

A typical color-invariant physical state is a linear combination of states 
of the form 



K 



N-'^l'al\-{k,)al^^{k,)...al:-{Km. 



(1) 



Here we sum over all possible values of the row and column indices. The 
superscript K is the integer sequence fci, A;2, • • • , ^c- The factor of in this 
equation serves as a normalization factor in the large- A^ limit. If we split K 
into two sequences Ki and K2 such that K1K2 = K, it follows that 



(2) 



where e{Ki) = if the number of integers in Ki larger than A^ is even, and 
e{Ki) = 1 if it is odd. 

A color-invariant observable is a linear combination of operators of the 
form 



75 = A^-^'^-'^«i?(^iX 



^2 



<rOa)<::?(ja-i)---<Oi)- 

Again we sum over all possible values of the row and column indices. / is 
the integer sequence ii, i2, . . . , 



and J is the sequence ji, j2, 



3a- 



The 



factor of A^ is chosen to ensure that once we act this operator on a state of 
the form given by Eq.(|l|), we will get another normalizable state. 

It is known that as A^ — >■ 00, if there were no fermions, the action of 7^ 
on \E'^ would be such that a segment of any cyclic permutation of K which 
is identical to J would be replaced with J rB|. This is the key reason why 



there is a connection between large- A^ matrix models and spin chains. With 
the inclusion of fermions, the situation is similar except that we have to be 
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careful with the signs of the terms: 



KiK2=K KiK2=K 
+ Y (^_iyiKiMK2)+eiK3)]^K2^IK3Ki 

KiK2K3=K 
+ Y (^_iyiKiHK2)^K2^IKi 

KiK2=K 

JiJ2=J KiK2K3=K 

In Eq.(^, a summation hke J2kiK2=k means that we sum over all possible 
pairs of Ki and K2 such that the concatenated sequence is K. A typical delta 
function like yields 1 if the sequence K is the same J, and otherwise. 

A crucial observation at this point is that Eq.(^, an equation in the 
context of quantum matrix models in the large- limit, can be paraphrased 
in the context of a quantum spin chain containing both bosons and fermions 
as follows. Consider a spin chain with c sites satisfying the periodic boundary 
condition. Each site can be occupied either by a boson or a fermion. There 
are possible bosonic states, and Ap possible fermionic ones. Let Ap{l), 
Ap(2), . . . , Ap{Ab) be the annihilation operators of these bosonic states, and 
Ap{AB + 1), Ap{Ab + 2), . . . , Ap^As + Ap) be the annihilation operators 
of the fermionic ones. Attaching daggers as superscripts to these operators 
turns them to the corresponding creation operators. Then a typical state of 
a spin chain can be written as 

ij"" ^A{{k,)Al{k,)---Al{K)\0). 

Define the Hubbard operator [19| as follows: 

x;^^AlmpU). 

It is readily seen that the action of a product of Hubbard operators on i/j^ is 

('_l~)[«(*a)+<=Oa)][<:{fcp)+e(fcp+i)H he(fcp+a_2)]^fcp+ti-l _ 

V / Ja 

\[e(ia-i)+e(ja-i)][e(fep)+<:(A.p+i)H he(A;p+a_3)] rA:p+a-2 

^_'^y(i{i2)+e{j2)]e{kp) ^kp ^kik2...kp-iiii2...iakp+akp+a+i---kc 
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for p + a — 1 < c. 

We are now ready to identify supersymmetric matrix models with super- 
symmetric quantum spin chains. Consider the following weighted sum of all 
cyclic permutations of the states of a quantum spin chain in the following 
way: 

Kik2=K 

where K2 may or may not be empty but Ki has to be non-empty. This 
obeys Eq.(^), too. Then after some algebra, it can be shown that 

= (^_lY<ia)MJa)Mjl) + eU2) + -+eija-l)] . 

(^_l'jHia-l)+e(ja-l)Mji)+t{j2) + -+e{ja-2)] . . . _ ]_ ^ [':(*2 )+e(j2 )] ) . 
c 

p=l 

because both the left and right hand sides of it satisfies Eq.(^. As the 
Hamiltonian of a typical matrix model is of the form Y^j j hj'^j, where only 
a finite number of /i/'s are non-zero, Eq.(^) provide us a way to transcribe 
the matrix model into a quantum spin chain with fermions and/or bosons, if 
for each non-zero hj the numbers of integers in / and J are the same. 

Let us give some typical examples to illustrate the idea above. Consider 
the following supersymmetric matrix model with = Ap = 2: 

matrix _ / ,34 , „,34 „,43 „,43 , „,12 , ,21 ,,12 „,21 



i^HubbTrd = - i7?2 + lit - 7S - lit + 734 + 734 - 743 ' 743, 

- (73\^ + 7?3^ + lit + itl) + {ill + ill + lit + itl) 

Eq.(BD tells us that the corresponding spin chain model is 



-HHubbard - - 2^ 2^ [-^p ^p+1 ~ ^p 
p=l (7=3,4 
n 

E( v31 v42 I \^24 , \rlA -v^23 , ^32 

y-^p -^p+i + Ap^;^ + Ap Ap_^i + Ap Aj 



32 vAl 
p+1 

p=l 

_ 1^42 1^31 _ x-13 1^24 _ TL-23 1^14 _ yAl t^32 \ 
^p+1 ^p ^p+1 ^p ^p+1 ^p ^p+1 ) 

n 

E\ ^ ( v2cr TAO-2 tact2 v2cr \ 

1^ \^p ^p^\- ^p -^p+lj 
p=l o-=3,4 



22 
P ■ 
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If we identify the states 1, 2, 3, 4 to be the vacuum state, the state with one 
spin-up and one spin-down electrons, the state with a spin-up electron only, 
and the state with a spin-down electron only respectively, we will see that 
the Hamiltonian in Eq.(^) is nothing but that of the Hubbard model |]2D . 
More specifically, let Cp,| and Cp,| be the annihilation operators of spin-up 
and spin-down electrons at the p-th site respectively. Then we can make the 
following identifications: 



Cp,| ^ ^p(3); Cp^i ^ ^p(4); and Cp^^Cp,! ^ 
The Hamiltonian can now be rewritten as 



Ap(2). 



Hubbard 



where np| = Cp^Cpi and = c^^^Cpi are the number operators for spin-up 
and spin-down states at the z-th site. This model was shown by Lieb and 
Wu to be exactly integrable using the Bethe ansatz |21[, and is believed to 
describe many condensed matter phenomena like high- Tc superconductivity, 
fractional quantum Hall effect, superfluidity, ... etc. 

It is not necessary for = Ap for a matrix model to be integrable. 
Consider another celebrated fermionic spin chain model, the t-J model 
r9[. Its Hamiltonian is given by 



E E 



+ U^Hp^Upi. 
p=l 



tj 



— i E E (Cpo-Cp+l^o- + C^p+l^a^^pcr) P 



n r 



J " r 1 

+ 2 E - I ('^pT + ^Pi) (^P+i,T^P+i,J,) 

p=l 



In this equation. 



P = Y[{l-np^npi) 

p=i 

is the projection operator to the collective state in which no site has 2 elec- 
trons. Sp is the usual spin operator: 



Sp 
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cy _ i_ c _ ct 
^ ~ 2i ' 

Sp = Cp|Cpi-; and 
Sp = Cp|Cp|. 

If J -C t, then this model is equivalent to the Hubbard model in the large-f/ 
limit. However, this model is in general different from the Hubbard model. 
When the values of J and t are such that J = 2t, we say that the model is at 
the supersymmetric point, and it is integrable. The corresponding integrable 
matrix model in the large- limit is 

i^r-^ = -t {ii! + ifi + lit + iti) - J {lit + iti - lit - ill) ■ 

Note that there are 2 fermionic states (States 3 and 4) but only 1 bosonic 
state (State 1) in this integrable model. 

To date many exactly integrable generalizations of the Hubbard model 
have been found |Q. We can use the same transcription procedure to write 
down the corresponding integrable supersymmetric matrix models in the 
large- limit. 

Besides considering matrix models in which there are only matrix degrees 
of freedom, we can also consider models, the matrix-vector models, in which 
there are both matrix and vector degrees of freedom. Such models describe, 
say, mesons in quantum chromodynamics or both quarks and squarks in 
supersymmetric Yang-Mills theory. For the sake of simplicity, we will confine 
ourselves to models in which there are only one column vector and one row 
vector degree of freedom. Let x'^ ^ind Xf^y where again runs from 1 to A^, 
be annihilation operators such that typical commutation relations hold. In 
particular, they obey the following two non-trivial commutators: 



Sill; and 



Ami ) A 



(7) 



The following formalism will still hold if we replace the commutators among 
X^'s and x^'s with anti-commutators. A typical color-invariant state is a 
linear combination of the form 



s 



K ^ N~(^^'y'x^^^ai-^ih)al7ik2) ■ ■ ■ alT^iK)xljO). 
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An observable is again a linear combination of jj introduced above. Never- 
theless, we find it convenient to introduce the following two kinds of opera- 
tors: 



and 



<,njiX02)---<rHja)x^". 

Actually an Ij is a finite linear combination of 7j's, and so is an rj, though 



it is no need for us to bother about the details here. (See Ref. [p5| for a 
comprehensive discussion of this.) 

Matrix-vector models correspond to spin chain models satisfying open 
boundary conditions. The relation between a 7 and Hubbard operators is 
almost the same as Eq.(|^), except that the summation index p runs from 1 
to c — a + 1 this time. It can be easily seen that the corresponding Hubbard 
operators for an Ij and an Vj are 

and 

(^_iy<ia-l)+e{Ja^l)Mjl)+e{j2) + -+eiJa-2)] . . . ( _ ]^ ) [e( j2) + e{i2 ) . 

Vhjl Y'^'2^2 Viaja 

^n~a+l^n-a+2 ' ' ' 

respectively. 

Let us see how the transcription is put into practice. Consider a matrix- 
vector model with the following Hamiltonian: 

- {ill + ill + 7lf + 7fl) + (ill + ill + lit + itl 
-U-il + (7! + ll + 272) 
-PT (^3' + -3 + 11 + rt) - p, (ll + rl + ll + rl 



It turns out that the corresponding spin chain model is nothing but a Hub- 



bard model with an open boundary condition pq, 27, 28 



n— 1 n 



(nil + riii) - (rill + - pi (ni^ + n„|) 



i=l 



(We remark that this particular Hubbard model is taken from Ref . [p^j . ) 

We have seen in the above discussion how to determine the integrability 
of a supersymmetric quantum matrix model in the large- limit — rewrite 
the Hamiltonian in terms of Hubbard operators, and then test if the Bethe 
ansatz is applicable by checking if the scattering matrix of the equivalent spin 
chain model satisfies the Yang-Baxter equation. As for the case of bosonic 
spin chains, we remark that the operators 7, I and r defined above form a Lie 
superalgebra The symmetry of the simplest integrable matrix model, 
the Ising matrix model, can be described by the Onsager algebra, which is 
a subalgebra of the cyclix algebra |T^. As a result, all conserved quantities 
of the Ising model are elements of the cyclix algebra. Likewise, we speculate 
that the Lie superalgebra describe the symmetry of supersymmetric matrix 
models, and may even help us determine the integrability of them. 

The support in part provided by the U.S. Department of Energy under 
grant DE-FG02-91ER40685 to us is acknowledged. 
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